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ABSTRACT 

The measurement of shape parameters of sources in astronomical images is 
usually performed by assuming that the underlying noise is uncorrelated. Spatial 
noise correlation is however present in practice due to various observational effects 
and can affect source shape parameters. This effect is particularly important for 
measurements of weak gravitational lensing, for which the sought image distortions 
are typically of the order of only 1%. We compute the effect of correlated noise on 
two-dimensional gaussian fits in full generality. The noise properties are naturally 
quantified by the noise autocorrelation function (ACF), which is easily measured in 
practice. We compute the resulting bias on the mean, variance and covariance of the 
source parameters, and the induced correlation between the shapes of neighboring 
sources. We show that these biases are of second order in the inverse signal-to-noise 
ratio of the source, and could thus be overlooked if bright stars are used to monitor 
systematic distortions. Radio interferometric surveys are particularly prone to this 
effect because of the long-range pixel correlations produced by the Fourier inversion 
involved in their image construction. As a concrete application, we consider the search 
for weak lensing by large-scale structure with the FIRST radio survey. We measure the 
noise ACF for a FIRST coadded field, and compute the resulting ellipticity correlation 
function induced by the noise. In comparison with the weak-lensing signal expected 
in CDM models, the noise correlation effect is important on small angular scales, but 
is negligible for source separations greater than about 1'. We also discuss how noise 
correlation can affect weak-lensing studies with optical surveys. 
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1. Introduction 

The measurement of source morphologies in two-dimensional images is a fundamental problem 
in astronomy. The measurements of shape parameters are usually performed while assuming that 
the underlying noise is uncorrelated. However, spatial correlation of the noise is always present 
to some degree, and can significantly affect the derived parameters. In experimental situations, 
noise correlation can be produced by various effects such as convolution of background light with 
a beam (or point-spread function), interferometric imaging techniques, CCD readouts, etc. 

This effect is particularly important for measurements of weak gravitational lensing. Weak 
lensing provides a unique opportunity to measure the gravitational potential of massive structures 
along the line-of-sight (for reviews see Schneider et al. 1992; Narayan & Bartelmann 1996). This 
technique is now routinely used to map the potential of clusters of galaxies (see Fort & Mellier 
1994; Kaiser et al. 1994, for reviews). Detections of the more elusive effect of lensing by large-scale 
structure have been reported by Villumsen (1995) and Schneider et al. (1997) in small optical 
fields. The search for a strong detection of the effect on larger angular scales is currently being 
attempted with present and upcoming wide-field CCDs in the optical band (e.g. Stebbins et al. 
1995; Kaiser 1996; Bernardeau et al. 1997), and with the FIRST radio survey (Kamionkowski 
et al. 1998; Refregier et al. 1998). The main challenge comes from the fact that weak lensing 
induces image distortions of only about 1%, and thus requires high-precision measurements of 
source-shape parameters. Correlated noise can produce spurious image distortions and correlations 
in the shapes of neighboring sources, and therefore must be carefully accounted for in weak lensing 
studies. 

The effect of correlated noise is also particularly important for radio surveys performed with 
interferometric arrays. In such surveys, images are produced by Fourier inversion of the visibilities 
from a set of antenna pair. As a result of the incomplete visibility coverage, part of the noise on 
the image plane contains long range correlations which can extend across an entire field. The effect 
of these correlations is particularly relevant in the context of recent large radio surveys such as 
FIRST (Becker et al. 1995; White et al. 1996) and NVSS (Condon et al. 1997). While the present 
analysis was motivated by our attempt to measure weak lensing by large-scale structure with the 
FIRST radio survey, most of the following is general and can be applied to any wavelength and 
imaging technique. 

Condon (1997) computed the errors in gaussian fit parameters in the absence of noise 
correlation. He also gave a semi-quantitative treatment of the noise correlation which relies 
on simulations. In the present paper, we generalize his approach to include a general analytic 
treatment of the noise correlation. We show how the spatial correlation of the noise can be 
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naturally characterized by the noise Auto-Correlation Function (ACF). We explicitly derive 
the effect of noise correlation on source parameters for a general two-dimensional fit and for a 
two-dimensional gaussian fit. We also consider correlations induced in the parameters of nearby 
sources. In particular, we focus on the ellipticity correlation functions, which are used in searches 
for weak lensing by large-scale structure. We apply our formalism to the case of the FIRST radio 
survey. After measuring the noise ACF for one FIRST coadded field, we compare the noise-induced 
ellipticity correlations to those expected for weak lensing by large-scale structure. We show that, 
while they are important on small angular scales, noise correlation effects are negligible for source 
separations greater than about 1'. 

This paper is organized as follows. In jg, we define the noise ACF and show how it can be 
measured in practice. In §|3j, we consider a general two-dimensional least-square fit. We derive the 
bias produced by the noise correlation on the mean, variance and covariance of the fit parameters, 
and the correlation of the parameters for two neighboring sources. In §||, we apply these results to 
the case of two-dimensional gaussian fits and compute the error matrix. We explicitly derive the 
ellipticity correlation function induced by the noise correlation. In §||, we consider the concrete 
case of the FIRST radio survey. We measure the noise ACF for a FIRST field, compute the 
induced ellipticity correlation function, and compare the latter to that expected for weak lensing. 
In §||, we discuss the implications of this effect for optical surveys and related searches for weak 
lensing by large-scale structure. Finally, §0 summarizes our conclusions. 



2. Characterization of the Noise 

2.1. Noise Auto-Correlation Function 

Let us consider a two-dimensional image with intensity /(x), where x = (xi,X2) is the pixel 
position. The total intensity can generally be decomposed into 

7(x) = S(x)+iV(x), (1) 

where S'(x) is the intensity of detected sources and N(x) is that of the noise. We take the term 
"noise" to broadly refer to any intensity which is not associated with detected sources. Depending 
on the context, it can include not only instrumental noise, but also the background light produced 
by undetected sources. 

We assume that, for each pixel, iV(x) has mean and variance given by 

(N(x)) = and (2) 
(A(x) 2 ) = a%, (3) 

respectively. The brackets denote an ensemble average which, in practice, is well approximated by 
a sky average. If necessary, the first equation can be enforced be subtracting a constant term from 
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iV(x). The noise is generally correlated from pixel to pixel. This is quantified by the noise ACF 
which we define as 

n(x°») = (N(x a )N(x b )), (4) 

where x ab = x b — x a . 

In the case of uncorrelated noise, the noise ACF becomes 

. . I a%, if x = (uncorrelated) 

Vi*) = \ n - (5) 

I I), otherwise 

In the continuous limit, i.e. in the limit of small pixel separation h, this can be written as 

?](x) ~ cr^r/i 2 (5^(x), (uncorrelated) (6) 

where 5^ is the two-dimensional Dirac-delta function. 

If the noise is gaussian and homogeneous (i.e. statistically invariant under translations), 
its statistical properties are completely characterized by f?(x). Generally, however, the noise is 
non-gaussian, and higher correlation functions are required for a full description. Nevertheless, 
since our analysis only involves the two-point correlation function r/(x), it is also valid in the 
non-gaussian case. In addition, the noise is generally not homogeneous. For instance, the 
sensitivity and beam shape can vary accross the field-of-view and thus cause the noise properties 
to depend on position. In this analysis, we assume that this effect is small and that ??(x) provides 
an adequate characterization of the noise, at least in an average sense. 



2.2. Practical measurement of the noise auto-correlation function 

The noise ACF, r/(x) is easily measured in practice. To do so, discrete sources are (iteratively) 
removed above a given threshold. Then, ??(x) is measured by computing the average of equation (Q) 
over pairs of pixels separated by r?(x). 

Figure |l| shows a portion of a field in the FIRST radio survey. The contrast was enhanced to 
make the noise more apparent. Stripe-like patterns in the noise are clearly visible and indicate 
that the noise is correlated. Figure |2| shows our measurement of the noise ACF for this field. 
Sources were iteratively removed by excising pixels with intensities above 4cr, and by also removing 
an area 5.4" (1 beam width) in radius surrounding each pixel. The measurement of rj(x.) was 
then performed by randomly choosing pairs of pixels with separations equal to x. The long range 
correlations reminiscent of the VLA antenna pattern are apparent. In §||, we discuss the case of 
the FIRST survey in detail. 

It is easy to show that, in the absence of correlation, the standard deviation, cr[r/(x)], of the 
measured ACF is 

a[rj(x)] = a 2 N /yfN pairs (x), (7) 
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Fig. 1.— Noise patterns for the coadded field 07210+29486E in the FIRST survey. A 9' x 9' region 
in the 46.5' x 34.5' of field is shown. The contrast for the noise was enhanced by clipping pixels 
with intensities higher than 2 mJy beam -1 . The residual rms intensity is ctn ~ .15 mJy beam -1 , 
for pixels at least 5.4" away from pixels with intensities greater than 0.6 mJy beam -1 . Apart from 
bright sources, "stripes" in the noise are clearly apparent. 
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Fig. 2. — Normalized noise auto-correlation function measured for the coadded field 07210+29486E 
in the FIRST radio survey. For clarity, only the inner 2' x 2' arcmin region is shown. The correlation 
function was normalized to 1 at x = 0. For this field, r/(0) = cr^r ~ 2.08 x ICT 8 mJy 2 beam -2 . For 
this measurement, the SNR for 77(0) is 60. 
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where N pa i rs (x.) is the number of pixel pairs with separation x. This provides a measure of the 
signal-to-noise ratio (SNR) of the ACF measurement, namely 



SNRMx)] = r?(x)/a[77(x)] = 




In the case of figure ^, the number of pairs was chosen to be N pa i rs = 3600 for all x. As a result, 
the SNR for r/(0) is 60. This high value of the SNR was achieved owing to the large number of 
pixels in the FIRST field. 



2.3. Noise Correlation from Beam Convolution 

As Condon (1997) remarked, a natural example of noise correlation arises when part of 
the noise is convolved with a beam or Point-Spread Function (PSF). The total noise is then a 
sum of a convolved component (e.g. background light which gets smoothed by the PSF) and an 
unconvolved component (e.g. Poisson noise from the detector). For simplicity, we ignore the latter 
component and suppose that the totality of the noise is convolved as 

iV(x) = J d 2 x'B(x-x')N(x'), (9) 

where iV(x) is the intrinsic noise which is assumed to be uncorrelated. The convolution beam 
B(x) is assumed to be normalized as J d 2 xB(x) = 1. 

As is easy to show using equation (|j) and equation (||) for N, the resulting noise ACF is then 

r?(x) = o^(x), (10) 

with 

/?(x ab ) = h 2 J d 2 xB(x - x a )£(x - x b ), (11) 
where, again, x afe = x b — x a . The function /3(x) can be thought of as the beam ACF. 

Figure ^ shows the dirty beam for a typical FIRST grid pointing. (The dirty beam is the 
image of a point source placed at the center of the field prior to CLEANing; see discussion in §[|). 
The coadded field of figure [l| is a weighted sum of 13 similar grid pointings, and therefore does 
not have a well defined dirty beam. It is nevertheless instructive to compare the noise correlation 
function of figure ^ to this dirty beam. Qualitatively, the resemblance is striking and reveals that 
most of the noise is effectively convolved with the dirty beam. 



3. General Fit with Correlated Noise 



There exist several methods for measuring the shape parameters of sources. The first method 
consists in fitting a model (usually a gaussian) to the two-dimensional source profile. Another 



-8- 



0.0 



0.2 



0.4 



0.6 



0.8 



1.0 



c 

» i— i 

e 

o 
u 
cd 



X 



Dirty Beam 



o 



l - 



-2 



1 1 1 1 1 1 1 

- 


1 1 | 1 1 1 1 1 1 1 1 1 | 1 1 1 1 1 1 1 1 1 | 1 1 1 1 1 1 1 1 1 

\ / 

- 




X ' 

: 




■ 









- 

III 1 1 



1 

X! (arcmin) 



Fig. 3. — Typical dirty beam for a grid pointing in the FIRST survey. This particular beam was 
derived using AIPS for the grid pointing 07000+50218 and was normalized to 1 at 6 = 0. This 
pattern is typical of VLA snapshot images in which limited coverage of the u-v plane is reflected 
in relatively intense sidelobes along the projected arms of the VLA. 
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possibility is to measure the multipole moments of the source. The latter method is usually 
implemented using various window functions to ensure proper convergence (e.g. Schneider & Seitz 
1995; Kaiser et al. 1995). The first method is more common in radio interferometric imaging 
(e.g. AIPS software package; see http://www.cv.nrao.edu/aips/), while the second is usually 
used with optical images (e.g. FOCAS, SExtractor packages; see Jarvis & Tyson 1981; Bertin & 
Arnouts 1996, respectively) 

Here, we only consider the effect of correlated noise on model fitting. An extension of our 



results to the multipole moment method is, however, straightforward. In § |3.1| , we derive the 
parameter corrections from the noise correlation for the case of a general two-dimensional x 2 -fit. 



We include corrections up to second order in the inverse SNR of the source. In §^2, we then 
compute the resulting bias in the mean (§ |3.2|) , variance, and covariance of the parameters (§ |3.3| 
In §|3.4| we compute the induced correlation in the parameters of pairs of neighboring sources. 



Finally, in §3.5, we consider general functions of the parameters. This analysis builds on the 



results of Condon (1997) to include an analytical treatment of the noise correlation. 



3.1. Least Square Fit 

Let us consider a sector of an image consisting of a single source superimposed on noise. 
We then consider a fit to the two-dimensional source profile with a function F(x;a), where 
a = (ai, ci2, •••) is a vector of parameters. In §||, we will consider the case where -F(x; a) is a two 
dimensional gaussian with a consisting of the gaussian normalization, peak location, and size 
parameters (see Eq. flolf). 

We assume that the fit is "good", i.e., that the source profile 5(x) is well described by 

S(x) ~F(x;a), (12) 

where a is the parameter vector that would be measured in an ideal measurement without noise. 
The total image intensity (Eq. jl]) thus becomes 

J(x) = F(x;a) + N(x). (13) 



We assume that the fit is performed using the method of least squares, without accounting 
for the noise correlation. This is usually the case for most fitting routines (e.g. AIPS). In other 
words, the fit is performed by computing the usual functional 

xVl-E ^-f™ 1 ' . d4) 

p a N 

where the sum runs over all the pixels in the image and x p are the position of the pixel centers. 
The best fit parameters a are found by minimizing x 2 ( a )- This minimum occurs when |^-(a) = 0, 
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that is, for the value of a which satisfies 



[F(x p ; a) - F(^; a) + N(^)} ^-(x?; a) = 0. 



(15) 



In the noiseless limit (JV — > 0), an obvious, if expected, solution is a = a. Note that in the 
most general case, there could exist more solutions if the system of equation is degenerate. We 
will not consider this complication here and assume that this solution is unique. 

In the presence of noise, the solution a will deviate from a. In most cases, the integrated 
source flux is much larger than the integrated flux of the noise. This is quantified by the source 
signal-to- noise ratio, SNR S which we take to be much larger than 1. We can thus treat the noise 
N(x p ) as a perturbation in equation (|T5|). For this purpose, we rewrite iV(x) as aiV(x), where a 
is a dimensionless parameter of the order of a ~ SNR^ 1 - We then expand the parameter vector in 
powers of a as a = a + aa' 1 ' + a 2 a^> + 0(a 3 ). 

After inserting these expansions in equation (EH), Taylor expanding, collecting terms in 
powers of a, and setting a = 1, we obtain 



a, = hi + DuPi + C imn P m P n + Eu km P m Qi k + 0(SNR^ 



where 



D 



F 
da 



Pijk 

Cimn 
Filkm 



^dF d 2 F , , 



da j dak 



F>k,rj ^^r,kj 



dF 

da. 



d F 1 



and 



Q 



E 



d 2 F 



daidaj 



( X P;a)AT(xP) 



(16) 

(17) 
(18) 

(19) 

(20) 
(21) 
(22) 

(23) 



Unless otherwise specified, the summation convention is assumed throughout this paper. For 
reasons which will become clear below, the matrix D is called the error matrix (see Condon 1997). 
In the oversampling regime, i.e. when the pixel spacing h is small compared to the source size, the 
sums above can be turned into integrals using the substitution 



^ - h 2 



2~P 



(24) 
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3.2. Bias in Fit Parameters 

We first focus on the bias induced on the source parameters by the noise correlation. By 
taking an ensemble average of equation (|l^), we find the mean of the parameter ai to be 

( ai ) = a, + D u (P t ) + C imn (P m P n ) + Eukm {PmQik) + 0(SNR7 3 ). (25) 

But from equation (|2|), we see that (Pi) = so that the first order term vanishes. Using the 
definition of t?(x) (Eq. [|j) we thus find 

( ai ) = ai + C imn V mn {0) + E akm U m i k {0) + 0(SNR7 3 ), (26) 

where 

^i(O) = EE^x^a^x^a)^), and (27) 

p q J 

The operand "(0)" was added for future convenience. In the matrices V and U, the noise ACF 
r/(x) is averaged over the pixel pairs weighted by parameter derivatives of the source profile. The 
noise correlation therefore produces a bias in the fit parameters of the order of SNR^ 2 . 

Note however, that this bias does not entirely disappear in the event of uncorrelated noise. 
Inserting equation (||) in the previous equations, we find that, for uncorrelated noise 

Vijip) = cr%Hij (uncorrelated), (29) 
hjk(0) = cr%E>ij 

Thus, the ensemble averaged parameter becomes 



Uijk(0) = cr%Bijk (uncorrelated). (30) 



(a { ) = ai - ^a 2 N B lkj DnD kj + 0(SNR s 3 ) (uncorrelated). (31) 

Therefore, even for uncorrelated noise, the x^-Rt parameters are biased to second order in SNR S . 
This is not surprising: while the x 2-n t is always unbiased if F(x; a) is linear in a (see e.g. Lupton 
1993), no such guarantee exist if it is non-linear. It is easy to check that (ai) is indeed unbiased in 
the linear case. However, most two-dimensional source models, including the gaussian model, are 
non- linear in their parameters. The correlation in the noise then makes this bias generally more 
pronounced. 



3.3. Parameter Variance and Covariance 



We can also compute the covariance of the fit parameters 

cov[cii,aj] = ((ai - (ai))(aj - (a,))) 



(32) 
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Inserting equation ([!(]) in this definition, we obtain 

cov[ ai , aj ] = D u D jk Vi k (0) + 0(SNR fi - 3 ). (33) 

In particular, the diagonal elements yield the variance (or squared error) of each of the fit 
parameters taken separately 

a 2 [ ai ] = D a D ik V lk (0) + 0(SNR s - 3 ), (34) 

where the index i is not to be summed over on the right hand side. 
In the case of uncorrelated noise, the covariance reduces to 

cov[aj,aj] = cr%Dij + 0(SNR~ 3 ) (uncorrelated), (35) 

in agreement with Condon (1997), and thereby justifying the term "error matrix" applied to D. 
The variance obviously becomes 

a 2 { ai ] = cr^Ai + 0(SNR7 3 ) (uncorrelated), (36) 

where, again, the index i is not to be summed over. 

3.4. Parameter Correlation Function 

Let us now consider two distinct sources, S 1 and <S 2 , in a single field. We assume that the 
sources are sufficiently distant from each other that the intensity of one source is negligible at the 
position of the other. If the noise in the field is uncorrelated, the parameter fit to each of the 
sources would thus be independent. But if the noise is spatially correlated, the parameters of the 
two sources will be correlated. This effect is particularly relevant for weak lensing studies in which 
one searches for pairwise correlations in the ellipticities of sources. 

To quantify the effect, we define the parameter correlation function as 

u,-(x 12 ) = ({a} - (a}))(a 2 - (a 2 ))) , (37) 

where a r is the fit-parameter vector for source S r , and x 12 is the separation between the (noiseless) 
source centroids. 



Inserting equation (flq) into this definition yields 



w 



« a (x 12 ) = DlD%V™{^) + 0(SNR s " 3 ), (38) 



where, 

^•-E^:( xP ^ r )^-(x p ;^), (39) 



„ dai daj 
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and 

W) = EE^i» 1 feA a M 3tM )- ( 4 °) 



da i ' daj 

p q 1 J 

The fit-parameter correlation is therefore of order SNR -2 . 

In the absence of noise correlation (Eq. ||), V^ 2 (x 12 ) becomes 



As is the case for a gaussian fitting function, the parametric partial derivatives, are often 
spatially localized. In this case, the above sum vanishes, if the sources are sufficiently distant from 
each other. Thus, for |x| much larger than the source sizes, and for well-behaved fitting functions 

io""(x)~0 (uncorrelated). (42) 

As expected, the fit parameters of the two sources are then uncorrelated if the noise is itself 
spatially uncorrelated. 



3.5. Functions of the Parameters 

We can also compute the effect of correlated noise on functions of the parameters. Let us 



consider a set of such functions, ij(a). In §4.2, we will consider an application were the tj's are the 
2 components of the ellipticity of the source. In the presence of noise, the value of the functions 
are perturbed as (see Eq. fill ) 

dt- 

ti(a) = ti(k) + -±(k)D jk P k + 0(SNR7 2 ). (43) 
The covariance of the t^s is thus 

= ^ ifc ^^(0) + O(SNR s - 3 ). (44) 

where we have defined 

Kij - ^)D kj . (45) 

In the absence of noise correlation (Eq. p9[), this covariance matrix becomes 

dt- dt- 

cov[t i ,t j }=a 2 N -^(a)-±( y k)D kl + 0( y SNRj 3 ), (uncorrelated) (46) 
oak oai 



in accordance with Condon (1997). 
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For a pair of two distinct sources, S 1 and S 2 , the correlation function for the tj's is easily 
found to be 



u^(x 12 ) ee ((hi* 1 ) - (t^ 1 ))) (t^s?) - ( tj ( a 2 



.1 

where 



K} k KlV^{± 12 ) + 0(SNR s - 3 ). (47) 



Kl^^Dl, (48) 

In the absence of noise correlation, and for distinct sources with well-behaved fitting functions 
(Eq. pH), this correlation function is simply 

wfj (x) ~ (uncorrelated) (49) 



4. Gaussian Fit 

In the previous section, we derived the bias, covariance, and correlation function of parameters 
for a general fit in the presence of correlated noise. In this section, we apply these results to the 
case which is of most practical interest, namely that of a two-dimensional gaussian fit. 



4.1. Gaussian Parametrization 

As the fitting function F(x;a), we consider the the following parametrization of a 
two-dimensional elliptical gaussian 

G(x, a) = ^ e -|(x-x") T A(x-x«) ; (50) 

where A is the amplitude, x a is the centroid vector, A is a symmetric positive-definite 2x2 
matrix which defines the shape and orientation of the gaussian, and the superscript T denotes the 
transpose operation. We choose the parameter vector to be 

a = (xl, x%, A, A n , A 12 , A 22 ). (51) 

Note that this parametrization is different from that of Condon (1997), and was to chosen because 
it is easier to relate to the parameters used in weak-lensing measurements. 

It is convenient to compute the first few multipole moments of G. With the above 
definition, we find 

?(x,a) = ttA\A\-\ (52) 
d 2 xxiG(x,ai) = xf, (53) 
d 2 x( Xi - xl){xj - x?)G(x,a) = (A-%, (54) 



G (o) 


~= 


<$) 


1 


G(°) 

r (2) 

G(°) 


~ G(°) 
1 

~ G(°) 
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where \A\ is the determinant of A. Note that 

j = a 1 = ur 1 



A 



22 



-A 



12 



v -A 12 A n J 

is exactly equal to the normalized quadrupole moments of G. It can be diagonalized as 

J = R(-a)( a Q 2 6 ° 2 )r(-«) T 

where a, b are the (la) major and minor axes, and a is the position angle measured 
counter-clockwise from the positive x-axis. The rotation matrix R is defined as 



(55) 



(56) 



cos (p sin (p 
— sin ip cos cp 



Inverting these relations yields 

(a 2 ,6 2 ) 
tan 2a 



Jn + J22 ± V( J n - J 2 2 ) 2 + 4J 2 2 

2Jl2 



<Al — ^22 



(57) 

(58) 
(59) 



4.2. Source and Noise Matrices 



We are now in a position to compute the source and noise matrices in equations (17-23). For 
this purpose, we first compute the partial derivatives of G. We find, 



dG 
da 



G ([(xt - xl)A xl + (x 2 - x%)A 12 ], [(x 2 - x%)A 22 + ( Xl - x$)A 12 ],A , 
-\{xi - xl) 2 , -(xt - xt){x 2 - x a 2 ), -i(x 2 - x'if 



(60) 



In the continuous limit (Eq. |24j), the components of the matrix H (Eq. |18|| ) have a closed 
form, i.e., 



H ~ 



ttA 2 
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where A 2 = (A11A22 + 2-Af 2 )- Note that the amplitude and shape parameters mix with each other, 
but not with the position parameters. 

By inverting the above matrix, we find the error matrix D to be 

/ 



D 



2h 2 ^\A\ 
irA 2 
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\A\- X A 22 


-\A\~ 1 A 12 














-\A\~ x A l2 


\A\~ l A n 
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AA 12 


AA22 
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4A 2 n 


4A n Ai 2 
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4A n A 12 


2(A n A 22 + A12 2 ) 


4A 12 A 2 2 








AA22 


AA 2 l2 


4A22A12 


^22 



J 



(62) 



The B matrix of equation (|lG|) can be computed in a similar way. The calculation for this 6 3 
component matrix is however cumbersome. Since it does not enter in the parameter correlation 
function on which we will now focus, we will not consider it further. 



The noise matrices P (Eq. p2[ ) and Q (Eq. ]23| ) can not be computed without knowledge of 
the noise ACF, rj(x). In practice, rj(x) has a complicated x-dependence. As a result, P and Q 
must be computed numerically using the direct measurement of ?](x), along with the expressions 
for the derivatives of G (Eq. In the next section, we study the effect of the noise correlation 

on the source ellipticity, a combination of the source parameters which is particularly relevant for 
weak lensing measurements. (By source ellipticity we mean the ellipticity of the observed source 
image and not the ellipticity in the source plane, as is common in the lensing nomenclature). 



4.3. Source Ellipticities 

The shear produced by weak gravitational lensing produces distortions in the shapes of 
background sources. The most direct way to detect this effect is to look for correlations in the 
ellipticities of sources. Measurements of the ellipticity involve integrals that are more spatially 
extended than that of the flux and position of a source. Consequently, we expect the ellipticity to 
be particularly sensitive to the spatial correlation of the noise. Here, we compute the magnitude 
of this effect on the ellipticities and on the ellipticity correlation function. 

A number of ellipticity measures have been used in weak lensing studies (e.g. Schneider 1995; 
Bonnet & Mellier 1995). Here, we will use the ellipticity measure considered by Miralda-Escude 
(1991), Kaiser & Squires (1993) and others, which has the advantage of being easily related to 



the gaussian parameters of equation (50). In this definition, the ellipticity is defined as a two 



component "vector" e whose components are given by 

^ _ - ^ _ 2jj2 ,gg. 

Jn + Jn + J22 

This can be rewritten in terms of the major axis, minor axis, and position angle of the source as 
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2 1,2 

e = fr^r (cos 2a, sin 2a). The component e±, sometimes written as e+, describes stretches along 
the x and y-axes, while e 2 , sometimes written as e x , describes stretches at ±45° with respect 
to these axes. In a coordinate system (x' 1: x 2 ) rotated by tp counter-clockwise from the positive 
xi-axis, the components of e are 

ej = R ii &p)e j , (64) 

where R is the rotation matrix defined in equation (|57|). This shows that e is not a true vector 
since it has a period of ir in <p. (In fact, e can be written as a symmetric traceless 2x2 tensor). 



The ellipticity can be rewritten in terms of the gaussian fit parameters (see Eq. |55|) as 

A22 - An -2A X2 fati , 

ei -^7T^' e2 = A^TM 2 (65) 

The partial derivative matrix of e is thus 

^ = 2(All+j422) ( A 12 -(A 11 + A 22 ) A 12 )- (66) 

Note that the position parameters, xf and x 2 , have decoupled from the ellipticity matrices and 
can thus be dropped. There are kept here for consistency. It is then straightforward to compute 
the K matrix (Eq. @). We find 

K~ 16h2 ^ ( A U -,4 22 \ 

TrA 2 (A u + A 22 )z\ A 12 \(An+A 22 ) A l2 J " 1 ' 

With this result and a measurement of the noise ACF r/(x), the covariance cov[ei,€j] and 
correlation function w?|(x) of the ellipticity are readily calculated using equations ([i4|) and (47). 

We now focus on «^j(x) which is of interest for measurements of weak lensing by large-scale 
structure. To be explicit, this correlation function is defined as 

™«(x 12 ) ee (e\e 

In practice, it is usually sufficient to consider the correlation of two sources which are identical 
and circular. In this simplified case, the unperturbed parameters for source S 1 and S 2 are 
a 1 = (A,x\,x 2 ,a~ 2 ,0,a~ 2 ), and a 2 = (A, x 2 , x 2 , a~ 2 , 0, a -2 ), where a is the (la) radius of the 
sources, A is their amplitude, and x\ gives their (unperturbed) centroid position. In this case, the 
ellipticity correlation function reduces to 

<M 2 )^^X2l\i /^x 1 d 2 x 2 X tJ (x 1 -x\ X 2 -x 2 ) e - (xl " Al) ^ x2 " 2) % ( x 12 ) + 0(SNR; 3 ), (69) 



where 

Tctfi & = ( mr - (mm? - ai?] zemm) 2 - ai?^ , 
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It is also convenient to consider the ellipticities measured with respect to axes parallel and 
perpendicular to x 12 , the vector connecting two sources ( Kamionkowski et al. 1998). Let us write 
this vector in polar coordinates as x 12 = (# 12 ,(^ 12 ), where 6 12 is the norm of the vector and tp 12 
the angle it substands with the xi-axis. From equation (p4|), the components of the ellipticities 
rotated into this coordinate system are e\ = Rij(2ip 12 )ej. The correlation function of the rotated 
ellipticities is then 

w%(9 12 ,<p 12 ) (efef ) = R^ 12 )R 3l (2^ 12 )wtK0 12 ^ 12 )- (71) 

Following the notation of Kaiser (1992), we identify Wj'J(x) = Ci(x), and ^^(x) = C2(x). We 
can also define a third independent correlation function u^2*( x ) — u; 2i( x ) = C3( x )> which should 
vanish if the noise is invariant under parity conservation. 

These correlation functions can then be averaged over ip 12 by defining 

= - f 'dpw$(e,<p). (72) 

J 7T Jo 

In the above notation, we simply write the azimuthally averaged correlation functions as 
w\\{0) = C^O), u)g(0) = C 2 {9) and w\l(9) = w r £{9) = C 3 {6). In the context of weak lensing 
studies, C\{6) and 62(0) can be directly related to the power spectrum of density perturbations 
along the line of sight (see e.g. Kaiser 1992). 

Figure |4j shows the rotated correlation functions for the FIRST coadded field of figures I and 
A source size of a = 2.45" and signal-to-noise ratio of SNR S = 5 were chosen. Figures ||, ||, and 
|?] show the corresponding azimuthally averaged correlation functions. In the next section, we will 
describe these results in detail and discuss their implications for weak lensing measurements. 



5. Application to the FIRST Radio Survey 

In this section, we apply the general formalism developed above to the specific case of the 
FIRST radio survey (Becker et al. 1995; White et al. 1996). The current version of the FIRST 
survey contains about 4 x 10 5 sources, and covers about 4.3 x 10 3 square degrees. The survey 
fields were observed at 1.4 GHz with the VLA in the B configuration. The 5a detection limit of 
the survey is about 0.75 mJy. The pixel size for the survey maps is 1.8", while the restoring beam 
is 5.4" (FWHM). The restoring beam is thus reasonably well sampled. This allows us to take the 
continuous limit (Eq. [^4j]) in our calculation. In addition, the relatively high angular resolution 
allows an accurate measurement of source morphologies. We are in the process of attempting to 
detect the effect of weak lensing by large-scale structure with this unique database (Kamionkowski 
et al. 1998; Refregier et al. 1998). 

Because of the observing time limitations for a survey of this magnitude, the observations had 
to be performed in the snapshot mode with integration times of 5 seconds. The interferometric 
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Fig. 4. — Ellipticity correlation functions for the coadded field 07210+29486E. Each correlation 
function C\, C2 and C3 is plotted as a function of the separation vector x in panels a, b, and c, 
respectively. For each source pair member, source signal-to-noise ratio was set to SNR S = 5, while 
the 1<t source size was set to a = 2.45". This corresponds to FWHM convolved and deconvolved 
diameters of 5.76" and 2", respectively. The clean beam for the FIRST survey has a FWHM 
diameter of 5.4". 
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Fig. 4b. — [See caption above] 



- 21 - 




Fig. 4c. — [See caption above] 
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data (or "UV" data) for each snapshot pointing was then Fourier transformed to construct a 
"dirty" image. The resulting pointing images were CLEANed using algorithms adapted from the 
standard AIPS deconvolution package. The final survey fields were produced by a weighted sum of 
the grid pointing images. At a given point, approximately 12 pointings contribute to the coadded 
field. Sources were then detected and their shape parameters measured, using a two-dimensional 
gaussian fit. 

Even though the grid pointing strategy was designed to optimize the signal-to-noise ratio 
of the coadded images, the final maps still suffer from the limited UV-coverage of the snapshot 
pointings. As a result, typical coadded maps contain small but noticeable "stripes" in the noise 
(see figure [j]). In other words, the noise is correlated. As we discussed above, this affects the 
source parameters derived by the fitting routine. In particular, the effect of correlated noise is a 
potential systematic effect in a search for correlations in the source ellipticities produced by weak 
lensing. 

5.1. Noise Auto-Correlation Function 

As a specific case, we considered the coadded field 07210+29486E in the FIRST survey. This 
field covers a solid angle of 46.5' x 34.5', with a pixel size of 1.8". Figure |l] shows a portion of 
this field with the contrast enhanced to make the noise more apparent. Bright radio sources 
are apparent in the right hand corner. In addition, the stripe-like nature of the noise is visible. 
For this analysis, pixels with intensities above .6 mJy beam -1 (corresponding to about 4cr) were 
excised from the field. The excised regions were then padded by 5.4" to avoid contamination of 
the noise by the source side wings. The residual noise has a standard deviation of &n — -15 mJy 
beam -1 . (Intensities are quoted in mJy beam -1 , were the beam refers to the restoring beam with 
a FWHM diameter of 5.4".) The median noise standard deviation for the FIRST survey is 0.14 
mJy beam -1 , with 95% of the FIRST area having a noise standard deviation less than 0.17 mJy 
beam -1 (White et al. 1996). This makes the coadded field 07210+29486E somewhat noisy, yet 
representative of the FIRST survey. 

Figure || shows our measurement of the noise ACF, r}(x), for this field. For this measurement, 
the SNR for 77(0) is 60 (see Eq. ||]). For convenience, r?(x) was normalized to 1 at x = in this 
figure. In addition to a prominent central peak, the ACF shows marked radial structures. These 
structures are reminiscent of the VLA antenna pattern and characterize the "stripes" which are 
visible in the noise. 

Each grid pointing is characterized by a dirty beam. (It is the Fourier transform of the 
UV coverage and is equal to the dirty image produced by a point source at the phase tracking 
center). In general, the dirty beam shape varies from one grid pointing to another due to differing 
geometries, flagged antennas, etc. As a result, coadded fields, which are composed of several grid 
pointings, do not have a well defined dirty beam. It is nevertheless instructive to qualitatively 
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compare our measurement of the noise ACF with a typical grid pointing dirty beam. Figure || 
shows such a dirty beam for the grid pointing 07000+50218. The qualitative resemblance with 
rj(x) is striking. This shows that a large fraction of the noise is effectively convolved with the dirty 
beam. 



5.2. Ellipticity Correlation Functions 

Of practical interest for weak lensing measurements, is the effect of the correlated noise on 
the ellipticity correlation functions. Figure |3] shows the rotated correlation functions Ci(x), C*2(x) 
and Ca(x) (Eq 71]) for the coadded field of figure |l| A source size of a = 2.45" was chosen. 



This corresponds to a FWHM diameter of 5.76", that is to a deconvolved FWHM diameter of 
2", for a restoring beam of 5.4". The source signal-to-noise ratio was set to SNR S = 5. These 
parameters correspond to both the detection and resolution limit of the FIRST survey and can 
thus be considered as "worse-case" parameters for a weak lensing search. 

The symmetric patterns of r/(x) (figure |2|) are also clearly apparent on figure ||. However, the 
patterns are different for each of the correlation functions. The first correlation function, Ci(x) 
exhibits a pronounced star-like pattern with radial elongation on each arm. On the other hand, 
C2(x) is mostly characterized by a negative correlation close to the center. Finally, C*3(x) exhibits 
star-like patterns but of much smaller amplitude than that of Ci(x) and Ca(x). The last fact is 
expected for a non-parity violating noise pattern. 

Figure || shows the azimuthally averaged correlation functions, C\{9), 6*2(0) and C%{9) as a 
function of 9. These are simply the radial profiles of the rotated correlation functions of figure |j. 
They were derived by Monte-Carlo averaging over 4 x 10 5 separation vectors x picked randomly 
within a radius of 50'. The thick lines correspond to a source size of a = 2.45", while the thin lines 
corresponds to a = 3.40". In both cases, the source signal-to-noise ratio was set to SNR S = 5, as 
before. For 9 <^ 0.4', C\{9) is positive while C2{9) is negative and drops off faster. Note that both 
Ci(9) and C2{9) are significantly non-zero well beyond the source radius a. This is of course due 
to the long-range features in r/(x) (figure ||). At all angles, C%{9) remains smaller and oscillates 
around 0. As we remarked above, this is expected for non-parity violating noise. 

A comparison between the thick and thin lines in figure [5] shows how the ellipticity correlation 
functions depend on the source size. For larger source sizes, C\{9) and 62(0) keep the same 
qualitative behavior but have a lower maximum amplitude and drop off more slowly. This is 
expected since larger source sizes effectively smooth out the sharp features in the noise ACF 
(figure ||; see Eq. 



Figure ^ shows the long range behavior of the ellipticity correlation functions. For angles 
larger than 9 = 10', a mosaic of 4 fields contiguous with 07210+29486E was used to compute 
the noise ACF rj(x). The error bars shown for C\{9) are the error in the mean derived from the 
Monte-Carlo sampling described above. Since the separation plane was undersampled, these error 
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Fig. 5. — Azimuthally-averaged ellipticity correlation functions for the same coadded field. Each 
correlation functions, C\, C2, and C3, are plotted as a function of the separation radius 9. A source 
signal-to-noise ratio of SNR S = 5 was assumed. The correlation functions scale as oc SNR^ 2 . 
The correlation functions are plotted for a source size of a = 2.45" (thick lines) and of a = 3.40" 
(thin lines). 



-25- 



bars provide a good estimate of the uncertainty in the correlation functions. As is apparent on 
this figure, all three correlation functions are consistent with for 9^1'. 

5.3. Consequences for Weak Lensing Searches 

As was demonstrated above, correlated noise affects the ellipticity correlation functions in the 
FIRST survey. We now compare the amplitude of this systematic effect with the signal expected 
for weak lensing by large-scale structure. 

Figure [?] shows the noise-induced correlation functions, C\(6) and C%(9) on a log-log scale. The 
absolute value of the correlation functions are plotted, with positive and negative values indicated 
by squares and stars, respectively. For comparison, the (absolute value of the) correlation functions 
expected for weak lensing in CDM models are also displayed. The four models correspond to the 
COBE-normalized CDM models in table 2 and figure 4 in Kamionkowski et al. (1998). Models 
1-3 correspond to a standard, tilted, and lambda CDM model, respectively. Model 4 corresponds 
to a standard CDM model with a smaller Hubble constant than that of Model 1. In all models, 
only linear evolution of density perturbations were considered. 

As we noted above, noise correlation produces significant biases in C\{9) and C2{0) for 9 <^ 1'. 
In this angular range, the noise-induced correlation functions exceed those expected for weak 
lensing for all models. (The inclusion of nonlinear evolution of density perturbations increases the 
estimates of the weak lensing signal for 9 <^ 10' by only a factor of a few and therefore does not 
significantly alter this conclusion; see Jain & Seljak 1997.) 

For angles 9 ^ 1', we showed above that the noise-induced correlation functions are consistent 
with zero, within the Monte-Carlo uncertainties. This can be seen by the alternation of squares 
and stars in figure g. The noise-induced lines in this figure should thus be interpreted as la 
upper limits in this angular range. (A larger number of Monte-Carlo simulations could of course 
reduce this upper limit even further but is computationally cumbersome.) We notice that, for 
1' ^ 9 ^ 50', these upper limits are almost one order of magnitude smaller than model 3, the most 
pessimistic CDM model. 

As we showed above, the noise-induced correlation functions depend only moderately on the 
source size (see figure |5|) . Most of the sources used in our weak lensing searches have sizes close 
to the resolution threshold (a ~ 2.45", or, equivalently, FWHM deconvolved diameter of 2"), 
and hardly ever exceed a ~ 5" (Refregier et al. 1998). As a result, the noise-induced correlation 
functions in such a survey will be close to that shown on figure 0. 

The Monte-Carlo realizations also allow us to compute the rms standard deviation a[Ci(9)] of 
the ellipticity correlation functions. For 10' < 9 < 50', we find a[C\{9)] ~ cr[C 2 {9)} ~ 1.9 x 1CT 3 . 
This is almost two orders of magnitude below the standard deviation produced by the intrinsic 
ellipticities of radio sources (a[C x (9)} ~ <r[C 2 (0)] « 0.4 2 = 0.16; see Refregier et al. 1998). Thus, 
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Fig. 6. — Same as the previous figure but this time focusing on the large-angle behavior of the 
correlation functions. The correlation functions are only plotted for a source size of a = 2.45". The 
error bars are the la errors in the mean derived from the Monte-Carlo integration. For clarity, 
they were only shown for C\{9). The correlation functions are all consistent with for 9 ^ 1'. 
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Fig. 7. — Comparison with the weak lensing correlation functions expected in CDM models. The 
smooth curves show the ellipticity correlation functions for the four COBE- normalized CDM models 
of Kamionkowski et al. (1998). The top and bottom curves for each model correspond to C\ and 
C2, respectively. The broken lines correspond to the ellipticity correlation functions induced by 
the noise correlation, for sources with a = 2.45" and SNR S = 5. In all cases, the absolute value of 
the correlation function is shown. To preserve a sense of the noisy behavior of the noise-induced 
correlation functions, positive and negative values are indicated by squares and stars, respectively. 
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statistical fluctuations in the noise-induced correlation functions are negligible for 6 ^ 1'. 

The estimation of the noise-induced correlation functions for 9 ^ 50' requires a large mosaic of 
coadded fields and is therefore computationally cumbersome. However, we expect the correlation 
functions to continue to be consistent with at large angles. In addition, a search for weak lensing 
will comprise an average over the whole survey, as opposed to an average over a few coadded fields 
in the present study. The dirty beam, and therefore t/(x), varies from one coadded field to the 
next. As a result, we expect any long-range features in the noise ACF to average out, leading to 
even smaller ellipticity correlations. 

We conclude that, in the FIRST survey, noise-induced ellipticity correlations dominate over 
the expected weak lensing signal for 6 <; 0.5', but are negligible for 9 ;> 1'. 

6. Correlated Noise in Optical Images 

While the effect is likely to be less severe than for intereferometric radio surveys, noise 
correlation may also affect source shapes in optical images. Most weak lensing studies are 
performed in the optical band and are thus potentially sensitive to this effect. Noise correlation 
in optical images can be produced by several effects. First, noise correlation can arise from the 
convolution of background light with the PSF. Because field distortions, diffraction spikes, etc, 
can produce anisotropy of the PSF, this can induce ellipticity correlations. In addition, CCD 
readouts and charge-transfer efficiency can potentially produce linear features in the noise ACF. 
Finally, observing and image processing techniques such as drift-scanning, shifting-and-adding, 
and "drizzling" (Fruchter h Hook 1997) can also produce features in the ACF. 

Measurements of the noise ACF in CCD images have recently been performed by Van 
Waerbeke et al. (1997), in the context of a novel technique to measure weak lensing shear with 
the noise ACF. In one of their fields, they find, at the center of the noise ACF, evidence for a cross 
pattern which they tentatively attribute to charge transfer efficiency and/or the shift-and-add 
procedure. Another application of their technique by Schneider et al. (1997) also indicates the 
presence of similar instrumental effects at the center of the noise ACF for the three CCD fields 
they considered. In one of the fields, the central region of the noise ACF is clearly elongated along 
one of the CCD axes. In another context, Vogeley (1997) measured the noise ACF of the Hubble 
Deep Field (HDF) in a search for fluctuations in the Extragalactic Background Light. Preliminary 
tests indicate that the extended wings of the PSF of the Hubble Space Telescope (HST) produce 
as much as 10-40% of this noise ACF signal. These wings are probably due to internally scattered 
light. Other effects such as flat-field errors could also contribute to the ACF but are likely to be 
less important. 

The effect of noise correlation should not be overlooked in weak lensing searches with optical 
images. We have shown that the effect of correlated noise is of second order in the inverse source 
SNR. One must therefore be cautious in correcting for systematic shape distortions using bright 
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stars. While this technique, which is standard in weak lensing studies in the optical, ensures 
proper correction of the effect of the convolution of the PSF with the sources, it would miss the 
effect of noise correlation. Indeed, for stars with high SNR, the effect of the correlated noise is 
negligible, while it could be significant for fainter galaxies used in weak lensing searches. 

It would thus be instructive to measure the noise ACF in optical images with high SNR and 
to compute the ensuing ellipticity correlations. This is one of the potential systematic effect in 
future high precision searches of weak lensing by large-scale structure with optical surveys (see 
Kaiser 1996; Bernardeau et al. 1997). In particular, the drift-scanning technique involved in the 
processing of the Sloan Digital Sky Survey could produce noise correlations which would need to 
be corrected for in future searches of weak lensing with this database (Stebbins et al. 1995). The 
effect of noise correlation could also be relevant for the recent detection of galaxy-galaxy lensing in 
the HDF (Hudson et al. 1997), since this detection involves measurements of ellipticities at small 
pair separations (a few arcsec) with the under-sampled PSF of the HST. 

7. Conclusions 

We have studied the effect of noise correlation on the shape parameters in two-dimensional 
images. The noise correlation is conveniently described by the noise ACF which can easily be 
measured in practice. We derived the magnitude of the effect for a general two-dimensional 
least-square fit. The noise correlation can produce a bias, and affect the variance and the 
covariance of source parameters. In addition, it can produce systematic correlation in the 
parameters of pairs of sources. We find that the effect is of second order in the inverse SNR of the 
sources. 

We applied these general results to the case of most practical interest, namely a two- 
dimensional gaussian fit. We computed the relevant matrices explicitly. In addition, we explicitly 
derived the systematic bias produced by this effect on the ellipticity correlation function of source 
pairs. This is particularly relevant for weak lensing studies. 

As a concrete example, we studied the effect of correlated noise on the shape of sources in the 
FIRST radio survey. We measured the noise ACF and found long range features which extend 
beyond the central maximum. We computed the resulting systematic effect on the ellipticity 
correlation function. We find that, in the FIRST survey, noise-induced ellipticity correlations 
dominate over the expected weak lensing signal for 9 <^ 0.5', but are negligible for 9 ^ 1'. 

We discussed the consequences of noise correlation for optical surveys. In optical images, 
noise correlation can arise from various effects such as the convolution of background light with 
the PSF, CCD read outs, shift-and-add preprocessing, drift-scanning, etc. Because the effect is 
quadratic in the source SNR, the effect could be overlooked if systematic distortions are monitored 
solely with bright stars. The effect of noise correlation could thus be important for searches of 
galaxy-galaxy lensing and of weak lensing by large-scale structure in optical surveys. 
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